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Measurements of ellipticities of background galaxies are sensitive to the reduced shear, the cosmic 
shear divided by (1 — k) where k is the projected density field. We compute the difference between 
, shear and reduced shear both analytically and with simulations. The difference becomes more 

important on smaller scales, and will impact cosmological parameter estimation from upcoming 
' experiments. A simple recipe is presented to carry out the required correction. 
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I. INTRODUCTION 



One of the most fascinating aspects of general relativity, and the first triumph for the theory, is that gravitational 
potentials can act as lenses for light from distant sources. The presence of large-scale structure and its associated 
potentials along the line-of-sight to distant galaxies implies that the images of most distant galaxies are slightly sheared 
compared to their intrinsic shapes. This shearing effect encodes information about cosmological distances and the 
J> ■ evolution of large-scale structure. For this reason gravitational lensing of background galaxies by large scale structure 
(cosmic shear) offers an excellent way to study the distribution of matter in the universe [j], 0. B iM ■ Measurements 
of the cosmic shear are already enabling us to constrain the dark matter abundance and clustering amplitude among 
other parameters E3> 0] • I n the future, large surveys may well uncover properties of dark energy, such as its 
abundance and equation of state [jj, IToL ITTl , and of neutrinos [l3L Il4| . This program will be successful only if we 

can make very accurate theoretical p~d~io^ M 

As experiments begin to go deeper and cover more and more sky, theorists must make sure that predictions are 
accurate enough to extract cosmological information in an unbiased fashion. There is a quantitative way to phrase this 
directive: the systematic errors on cosmological parameters induced by theoretical uncertainties should be significantly 
I smaller than the anticipated statistical errors. Since the latter hover near the percent level for the most ambitious 
O | experiments, theorists clearly have their work cut out for them. 

-£3 ■ Here we consider one correction to the standard theoretical predictions, the effect of reduced shear [ilj ITtI Il8l | . The 
observed ellipticities of galaxies (with two components gi for each galaxy) are often used as estimates of the cosmic 
shear (ji), but in fact they are sensitive to the reduced shear: 

, M<i (i) 
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where k is the convergence or roughly the projected densityfield. We analyze the difference between shear and reduced 
shear both analytically and with numerical simulations |l9t l20j|. focusing on various two-point functions. Then we 
map out the region in experiment-space where the effects of reduced shear need to be included. Outside of this region, 
the canonical prediction - which neglects the 1 — k denominator - is sufficient. 

Throughout we assume a flat, ACDM cosmology. Our base model has f2 m — 0.28, flbh 2 — 0.024, h = 0.7, 
n = 1, and erg = 0.9. We will let the galaxy density and sky coverage of surveys vary, but we limit ourselves to all 
background sources at z = 1. Our quantitative results will change slightly for higher redshift sources, but our two 
major conclusions - that we can compute these corrections accurately and that we have to compute them if we want 
to extract cosmological parameters - are only strengthened (since the effect is larger) for higher redshift sources. 



II. SHEAR TWO-POINT FUNCTIONS 



Here we briefly review a variety of definitions relating to cosmic shear and its statistics. Cosmic shear can be 
represented by two numbers at any point in space, 71 and 72. Similarly, the ellipticity of a background galaxy can 
be described by g\ and g^. The latter are measurable, while the former are related in a straighforward way to the 
projected gravitational potential and therefore are simplest to compute given a cosmological theory. On average, all 
these components are zero; however their two-point functions contain information about cosmic fluctuations. 



2 



We focus here on four sets of two-point functions. 

• Smoothed Variance Defining 7 = 71 + 172, make a map of 7 smoothed over a square pixel of side 8. The 
variance of the smoothed shear field, 7, is then 

(\i 2 \m = (il)(e) + {ilm (2) 

with a similar definition for the galaxy ellipticities, g. Note that " ~ " denotes an average over a local square 
while angle brackets denote an average over the sky. 

• Aperture Mass At any sky position Oq, this is a 2D integral over the tangential shear, j t = — 7i cos(2</>) — 
72 sin(20), where 4> is the angle between Oq and a fixed x-axis. The wei ghti ng function in the integral depends 
on a smoothing scale 8. Here we use the smoothing function defined in [lfij . The average value of M ap is zero, 
but its variance as a function of smoothing scale contains information about the underlying fluctuations. 

• Correlation Function Since there are two components of shear, there are in principle 3 separate correlation 
functions (-fi(9 )-fj(8 +9)) averaged over all positions 8 . These depend on the angular difference 8 = \8\. Here 
we focus on the combination £(#) = (-f(6o)j* (0o+O)) . 

• Angular Power Spectrum Write the 7j field as a sum of coefficients times spherical harmonics. In the small 
angle limit in which we will work, this is equivalent to a Fourier transform, 7;(1). One linear combination 
of these Fourier coefficients (the so-called "B-mode" ) vanishes if the underlying fluctuations are due to scalar 
perturbations; the other, 

E(l) = £0-^(1)^(1), (3) 

is sensitive to the projected gravitational potential. Here is the 2D anti-symmetric tensor £12 = —£21 = 1; 
and the trigonometric weighting functions are 

Ti(l) = -sin(20 ; ) 

T 2 (l) = cos(20,) (4) 

where <pi is the angle of 1 with a fixed x-axis. The angular power spectrum is roughly the variance of these 
Fourier coefficents, 

(E(l)E(l')) = (2ir) 2 5 2 (l + l')C l . (5) 



Each of these two-point functions can be computed from a simulation using either shear or reduced shear. Thus, for 
example, we can measure in a simulation both and Cf and find the difference between the two. The simulations 
we use to compute these functions are described in Ref. lj|. We can also compute the two-point functions semi- 
analytically. The two-point functions of 7 can be ex pres sed in terms of integrals over the 3D matter power spectrum, 
which has been extremely well-studied |2lll2^l23ll24Ll25| . The two-point functions of g can be computed perturbatively 
by expanding the 1 — k denominator around k = 0. In §111, we write the shear two-point functions as integrals of the 
3D matter power spectrum (these expressions are well-known 0,0) and the reduced shear corrections in terms of the 
3D three-point function, the matter bispectrum. 



III. PERTURB ATIVE CALCULATION 



The simplest two point function to compute is the angular power spectrum. The Fourier transformed shear can be 
expressed in terms of the projected gravitational potential 

7,(1) = -e tf T i (1)^(1) (6) 

with 

m [°° ^w(x) I ^rHVx, h; x). (7) 

Jo X J 2tt 

Here \ is comoving distance; the lensing kernel W — 2%(1 — x/Xs)0(Xs — X) with Xs the distance to the source galaxy 
and O the Heaviside step function. 
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Inserting Eq. © into Eq. to get E(V), multiplying E(\) by E(V), taking the expectation value, and then 
integrating over 1' in the Limber approximation leads to 



I 4 f d X W 2 ( X ) p n/ , 

4 7 — x 5 — p * (//x;x) 



(8) 



where P$ is the 3D power spectrum of the gravitational potential. The top panel of Fig.^shows the power spectrum 
computed in this fashion as compared with that measured in simulations. Agreement is excellent, confirming earlier 
work [2?i 1271 . The one aberrant point on small scales in the simulations is close to the Nyquist frequency, so power 
in the simulation is artifically suppressed. 
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FIG. 1: Results from simulations (points with error bars) and the perturbative calculation (smooth curves) for the zero order 
power spectrum (top panel) and the fractional change due to reduced shear (bottom panel). 



Eq. ||SJ) is an expression for the power spectrum of cosmic shear, C ; 7 . To lowest order, when k — > in Eq. l|T]l. the 
power spectrum of the observable reduced shear is equal to this. We can perturbatively compute the correction to the 
reduced shear: to leading order, <?W = 7 and to second order — 7/t. Therefore, the correction to the two-point 
function due to reduced shear is 

5(E(l)E(l')) = wuTiWW)®™ (1)# 3 (1')> + (1 - 1') (9) 
Plugging in for E and using f-ijtjk = —dik, we have 

6(E(1)E(V)) = dL Ti (l) J 0^Ti(h)ll(l - li) 2 WliMl " hMl')) + (1 - 1') (10) 

Using Eqs. (19,20) and (22) in Ref. [23, we can reduce this to 

5Cl = ^P-J ^L Ti (h)ll(\ - h) 2 B%h, 1 - I X) -1) (II) 

where B K is the bispectrum of the convergence. Just as the power spectrum of the convergence can be written as 
an integral of the 3D power spectrum along the line-of-sight (Eq. ©), the 2D bispectrum is an integral of the 3D 
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bispectrum |29| : 

-Z 6 f°° W 3 (y) 

B K (h,l 2 ,h) = — / dx^p-B* (l l X , h/ X , l 3 X • (12) 
8 Jo X 

The 3D bispectrum, has been computed analytically on large scales and measured on a wide range of scales 
in simulations [24j. An accurate fit to the N-body results was introduced in Ref. [3(j; we use this fit to compute 
SCi, the difference between cosmic shear and reduced shear power. The power spectrum, P$, which is needed to 
compute B$, was computed using the publicly available Halofit code 2Sj which has also been calibrated by numerical 
simulations. Fig. ^ shows the results of this perturbative calculation and of a similar measurement from simulations. 
The perturbative results arc in excellent agreement with the simulations. This is extremely encouraging because it 
offers an easy way to include reduced shear corrections without resorting to expensive simulations. 

The conclusion that reduced shear differs from cosmic shear most significantly on small scales follows from per- 
turbation theory. On large scales, fluctuations in 7 and k are small; since the difference between cosmic shear and 
reduced shear is higher order in these perturbations (cx K7), it is very small on large scales. 

The angular two-point functions described in §11 can all be expressed as integrals over the power spectrum. The 
smoothed variance is 

d^DW = ^af ^joM/2).?o(M/2)Q. (13) 

This variance can be computed for either shear or reduced shear. The difference between the two is the same integral 
over SCi from Eq. (|ll|l . The other spatial functions are 

(O^^f^ (14) 

and 

m = ^- r dllJ o( w ) c i- ( 15 ) 

2tt Jo 

We have computed these three functions from simulations and perturbatively; the results are shown in Fig. [2J The 
simulations and perturbative calculations agree extremely well, as do the corrections. 



IV. IMPACT ON COSMOLOGICAL PARAMETERS 



When docs one need to include the effects of reduced shear when comparing models with observations? Neglecting 
these effects results in an incorrect prediction for the power spectrum; we computed a correction, SC), above. This 
incorrect prediction propagates to an incorrect estimate of the cosmological parameters, pi, or a bias. The bias on 
parameter i is 

b, ee p\™ - Pi = Fr/ wig^SQ. (16) 

Here wi is the weight, the inverse variance of the measurement, and F is the Fisher matrix 

dCidCi , , 



The variance depends on experimental specifications: sky coverage and depth/resolution. Specifically, 

^ = (q + ffi-V (is) 



(21 + l)/sky V n off 

The fraction of sky covered is / s k y , while the rms of the intrinsic ellipticity of galaxies, (7i„ t } 1 ^ 2 , is set to 0.25 [Toll3l|. 
and n e g is the effective galaxy density which depends mainly on the depth and resolution of the experiment. 

We can compute the bias induced by neglecting the difference between shear and cosmic shear for any experiment. 
For concreteness, we allow three cosmological parameters to vary: the normalization of the power spectrum, ag] the 
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FIG. 2: Angular two-point functions from simulations and perturbatively. 

0.022 

0.02 

0.018 

.5 0.016 
PQ 

0.014 
0.012 
0.01 




10 



100 



1000 



n eff (arcminute -2 ) 



FIG. 3: The bias due to neglecting reduced shear corrections as a function of galaxy density (per square arcminute for a 
1-component rms shear of 0.25) for the three cosmological parameters varied. The bias on as is negative. 
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shape of the power spectrum, T ~ £l m h; and the matter density, Vl m . Figure [3] shows the bias induced in these 
parameters by neglected the cosmic shear corrections. This bias was computed including data out to I = 3000. The 
cut-off is necessary because baryons affect the theoretical predictions on small scales [2M El , and it is very difficult 
to predict these effects accurately. 

Since the reduced shear corrections on scales I < 3000 are of order a few percent (Fig.^l, it is not surprising that 
the induced biases on the parameters are also of order a few percent. Note from Eq. Ijf 6(1 that the bias scales as 
F~ 1 wi. Since wi oc f~js., and F~ x scales as / s k y the bias is nearly independent of sky coverage. Therefore, the bias 
depends only on the galaxy density. At very large density, the weight wi becomes independent of galaxy density, since 
shape noise due to intrinsic ellipticity is inversely proportional to galaxy density. In this limit, cosmic variance - the 
first term in parentheses in Eq. I|I8|I - becomes the dominant source of noise. The largest bias is to the normalization 
parameter erg. 

How important is a one to two percent level bias in the cosmological parameters? We must compare the bias 
to the anticipated statistical error in an experiment. If there were only one free parameter, this comparison would 
be straightforward. With several parameters, we must ask which statistical error should be used as a baseline: the 
error on as for example which accounts for the uncertainty in all other parameters (this is called the marginalized 
error) or the error expected if all other parameters are held fixed? We argue that the latter should be used. To 
understand why, consider the case with two parameters. The expected I-er constraints trace out an ellipse in this 
plane. If the parameters are degenerate 1 , this ellipse is very elongated, and the expected constraints on either of the 
parameters individually will be very weak. If these marginalized errors are used, then the ratio (bias/error) will be 
quite small. A bias much smaller than the marginalized statistical error, however, can easily induce an analyst to 
estimate the parameter as lying outside the allowed ellipse. If we instead use unmarginalized errors as the baseline, 
then bias/statistical error ratios smaller than one are more likely to keep the parameters within the allowed ellipse. 
Fig. 21 shows the ensuing ratio of bias to statistical error for the most severely affected cosmological parameter (as) 
as a function of survey width (/sky) and galaxy density. 
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FIG. 4: The ratio of the bias on as due to neglecting reduced shear to the statistical error on as- The ratio depends on the width 
and galaxy density of the survey. Several current and future surveys are shown for orientiation: CTIO Lensing Survey 0; 
Supernova Acceleration Probe (SNAP) (34|; Dark Energy Survey (DES) http://www.darkenergysurvey.org ; Large Scale 
Synaptic Survey (LSST) |35J. If the ratio is larger than one, neglecting reduced shear corrections biases the parameters by an 
amount greater than the unmarginalized statistical errors. 



A smart analyst tries to choose non-degenerate parameters; in that case, the marginalized and unmarginalized errors are equal. 
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Alternatively, we can transform to a basis in parameter space where the Fisher matrix is diagonal and then rescale 
so that the errors on the parameters are unity. In this basis, magnitude of the bias is a simple measure of how far 
outside (or inside) the 1-sigma contours the bias leads. Again a bias greater than unity means the effect needs to be 
accounted for. Figure shows the bias in this rotated/rescaled basis. For example, neglecting the bias in the dark 
energy survey leads to an incorrect estimate of parameters at the 3-sigma level. 



Rotated/Rescaled Bias 




FIG. 5: The bias for rotated and rescaled parameters such that the anticipated errors are diagonal and unity. That is, the 
1-sigma contours in the 3D parameter space map out a sphere of radius one in this basis. The amplitude of the bias is shown 
here as a function of galaxy density and sky coverage. 

Both Figures 0] and El illustrate that the largest of current experiments are already entering a regime in which the 
effects of reduced shear must be considered and future experiments will certainly need to account for it. Relying 
on Fig. to justify neglecting the effect for current experiments is a little dangerous since it does not account for 
priors, which can reduce errors while leaving the bias unaffected. So the most conservative approach is to use Fig. 01 
as a guide to the impact of reduced shear. In that case, virtually all upcoming experiments will need to add in this 
correction when extracting information about cosmological parameters. 



V. OTHER CORRECTIONS 



Reducing the shear is not the only correction that needs to be applied to weak lensing spectra [id Il9t l'2("il| . For 
instance, Hu and Cooray |36j | have computed perturbative corrections to the power spectrum which account for 
the Born approximation and lensdens coupling. These corrections are an order of magnitude smaller than the one 
considered here. To understand why, recall that SCi from Eq. © comes from considering the product of a second 
order perturbation with a first order perturbation, x g^ 2 \ We have considered the corresponding terms for the 
beyond Born and lensdens corrections and found that they vanish. Therefore, the first non-vanishing corrections are 
the ones considered in Ref. [3(|: those of order g^g^ or g^g( 2 '. We have also not considered the effects of source 
clustering [33 . There is some indication j3^,[3!| that this may also induce percent level changes in the power spectrum 
on small scales. If so, these would need to be included as well. Accounting for this coupling in a simulation requires 
input from a real 3D galaxy catalogue. 
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VI. CONCLUSIONS 



Deflection of light rays by gravitational potentials along the line of sight introduces a mapping between the source 
and image plane. The Jacobian of this mapping defines the shear and convergence as a function of position on the 
sky. In the absence of size or magnification information neither the shear nor the convergence is observable, but only 
the combination g = 7/(1 — k), known as the reduced shear. On small scales, where lensing surveys get much of their 
constraining power, this must be taken into account when predicting the observables. 

We have studied the difference between cosmic shear and the reduced shear. Our main conlusions are: 

• The perturbative calculation of reduced shear, Eq. Ijllfl . agrees well with numerical simulations. This is not too 
surprising, since the lensing is weak and we use many components themselves fit to N-body simulations, but it 
gives us confidence that we can use these calculations in making predictions or fitting to data. 

• The effects of reduced shear are on the threshhold of becoming very relevant. As depicted in Figs. 0] and |5J 
upcoming surveys will need to account for reduced shear when extracting cosmological parameters. 

SD was supported by the DOE and by NASA grant NAG5-10842. MJW was supported in part by NASA and the 
NSF. CAS was supported in part by the Kavli Institute for Cosmological Physics at the University of Chicago and 
by grant NSF PHY-011442. The simulations were performed on the IBM-SP at NERSC. 
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